Macroscopic Klein Tunneling in spin-orbit coupled Bose-Einstein Condensates 
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We propose an experimental scheme to detect macroscopic Klein tunneling with spin-orbit cou- 
pled Bose-Einstein condensates (BECs). We show that a nonlinear Dirac equation with tunable 
parameters can be realized with such BECs. Through numerical calculations, we demonstrate that 
macroscopic Klein tunneling can be clearly detected under realistic conditions. Macroscopic quan- 
tum coherence in such relativistic tunneling is clarified and a BEC with a negative energy is shown 
to be able to transmit transparently through a wide Gaussian potential barrier. 
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I. INTRODUCTION 

Shortly after the relativistic equation of electron was 
established by Dirac, Klein used it to study electron scat- 
tering by a potential step and found that there exists a 
nonzero transmission probability even though the poten- 
tial height tends to infinity [1], in contrast to the scatter- 
ing of a non-relativistic particle. This phenomenon has 
been referred to as Klein tunneling (KT). KT is an intrin- 
sic relativistic effect and is interpreted as a fundamental 
property of Dirac equation that particle and antiparticle 
states are inherently linked together as two components 
of the same spinor wavefunction ,2,] . 

This unique scattering process has attracted lots of in- 
terest over the past eighty years but failed to be directly 
tested by elementary particles due to the requirements 
of currently unavailable electric field gradients [sj. In- 
terestingly, the dynamics of particles in some systems, 
such as electrons in graphene Q and trap ions 0,1^1 etc., 
may be described by effective relativistic wave equations 
and have been proposed to observe such relativistic tun- 
neling. Ultra-cold atoms in optical lattices @ and light- 
induced gauge fieldsH are also able to behave as rela- 
tivistic particles [a, l9|. Recent experiments in graphene 
hetero junctions [ICt llll ] have provided some indications 
for KT. However, the existence of disorders and interac- 
tions in these solid-state systems makes it hard to realize 
full ballistic scatterings. In addition, it seems hard to un- 
ambiguously observe KT in graphene since it is a typical 
two-dimensional (2D) system, while scattering in a 2D 
system is a combination of perfect transmission for nor- 
mally incident particles (a relativistic effect) and expo- 
nentially decay tunneling for obliquely incident particles 
(a non-relativistic effect). Moreover, KT as well as Zit- 
terbewegung effect have been experimentally simulated 
with the trapped ionsQ. 

In this paper we propose a feasible experimental 
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scheme to observe macroscopic KT with spin-orbit cou- 
pled BECs [n m. We demonstrate that a one- 
dimensional nonlinear Dirac equation (NLDE) with tun- 
able parameters can be realized with a spinor BEC in the 
presence of a light-induced gauge field. Through numer- 
ical simulations, we demonstrate that a macroscopic KT 
can be observed under realistic conditions. The simple 
configuration of gauge field, in combination with control- 
lable dimensions, interactions and potential barriers may 
provide us with a clean and tunable platform for investi- 
gation of interesting relativistic tunneling effects. 

We investigate the relativistic tunneling of a macro- 
scopic quantum object by comparing the transmission 
coefficients between a BEC in the absence of interactions 
and an incoherent ensemble average of non-condensed 
atoms. In addition, we find that a realistically weak 
interaction between atoms slightly affects the transmis- 
sion coefficients. The main feature of a BEC is that all 
atoms in the BEC are in the same state and in the same 
phase and thus the BEC can be considered a macroscopic 
object. So the tunneling of a BEC we study is the co- 
herent scattering of a macroscopic object. Tunneling in 
the former shows a distinct difference in relativistic ef- 
fects between macroscopic objects and the ensemble av- 
erage of some microscopic particles, while KT has been 
studied previously only within a single-particle scenario. 
We also present another unexpected result: that a BEC 
with a negative energy can almost completely transmit 
through a Gaussian barrier. Since KT is a relativistic 
phenomenon associated with an anti-particle in the po- 
tential, our proposed spin-orbit coupled BEC can mimic 
a macroscopic 'anti-BEC (a super-atom made from 'anti- 
atoms'), at least in a scattering problem. Therefore, the 
mimicked 'anti-BEC may open the possibility of explor- 
ing exotic relativistic effects of a macroscopic body (even 
for very large antimatter), in contrast to the conventional 
wisdom that relativistic effects are only clearer for a mi- 
croscopic particle. 

The paper is organized as follows. In Sec. II we pro- 
pose an approach to realize a spin-orbit coupled atomic 
gas through a A-level configuration, and then demon- 
strate that the dynamics of the atoms should be described 
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by the NLDE when the atoms are condensed into a BEC. 
In Sec. Ill we show that the region for KT of a single 
atom can be reached in experiments. Then we demon- 
strate in Sec. IV that KT of BECs can be clearly ob- 
served. We also clarify the macroscopic quantum coher- 
ence in such relativistic tunneling and show that a wide 
Gaussian potential barrier is transparent for a BEC with 
a negative energy. In Sec. V, we present our discussion 
and conclusion. In the Appendix, we briefly review the 
numerical method to calculate the transmission coeffi- 
cient of a single atom scattered by a Gaussian potential. 




FIG. 1: (Color online). Schematic illustration of the system, 
(a) Atom with a A-level configuration interacting with laser 
beams characterized by Rabi frequencies Qi, Q2 and a large 
detunning A. (b) Configuration of laser beams to realize a 
Dirac-like equation by the lasers ili ,Q2 and an effective Gaus- 
sian (square)-shaped potential induced by another laser beam. 
Atoms are confined in a ID waveguide along the y axis and 
scattered by the potential. 



II. REALIZATION OF A NONLINEAR DIRAC 
EQUATION WITH COLD ATOMS 

The Dirac equation with tunable parameters can be 
realized with ultracold atoms through two approaches 
[1, H 0- Similarly to graphene, it was proposed that 
low-energy quasiparticles in a honeycomb optical lattice 
should also be described by the relativistic Dirac equation 
0. On the other hand, the Hamiltonian of cold atoms 
(without optical lattices) with a certain spin-orbit cou- 
pling, which can be achieved with synthetic gauge fields, 
is a Dirac Hamiltonian when the wave number of the 
atoms is much smaller than the wave number of the laser 
beams. It is demonstrated that the required spin-orbit 
coupling can be realized though a tripod level configu- 
ration [1, Q ■ In this paper we proposed that a A-level 
configuration is also feasible for use in the realization of 
the Dirac equation. 

Let us consider the motion of bosonic atoms with mass 
m in the y-z plane, with each having a A-level struc- 
ture interacting with laser beams as shown in Fig. 1. 
The ground states |1) and |2) are coupled to an excited 
state 1 3) through laser beams characterized respectively 



and 



with the Rabi frequencies fli = J7 cos(Kj,2/)e 
n2 = rJsin(K2,?/)e*(^-'"-^\ where n = ^l^iP + l^2p. 
As shown in Fig.l (b), the Rabi frequencies ili and ft2 
can be realized, respectively, with a pair of lasers fii± = 
^nexp[i{—KzZ ± Kyu)] and fl2± = hi^exp{i[—KzZ ± 



(KyU + TT / 2)]} , where Ky = ncostp and = Ksimp with 
K being the wave number of the lasers and tp being the 
angle between the laser and the y axis. The Hamilto- 
nian of a single atom reads H = ^ -I- V{r) + Hj, where 
V{r) = E^=i(^T(r) + Vb{Y))\j){j\ denotes the full exter- 
nal potentials (including the trapping potentials Vt and 
the scattering potential Vb) and the interaction Hamilto- 
nian Hi = ;iA|3)(3| - {Y.l=i ^f^j|3)01 + h.c), with A 
as the detuning. Diagonalizing Hi yields the eigenval- 
ues a{[A- VA2 -|-4172]/2,0, [A-FVA2 + 4f72]/2}. In the 
large detuning case, the two eigenstates corresponding to 
the first two eigenvalues span a near-degenerate subspace, 
and can be considered a pseudo-spin with spin-orbit cou- 
pling induced by a gauge potential 0, Under this 
condition we obtain the effective Hamiltonian 



H 



pI 



2m 



Vy(7yPy + VzazPz+lzCrz+VT + Vb, (l) 



where 



and 7. = 



fi^/A^)^^] -I- In the derivation, we have dropped an 
irrelevant constant and assumed that the potentials V^(r) 
are spin- independent. Furthermore, the atomic gas can 
well be confined by a ID optical waveguide along the y 
axis [eI, so we may further restrict our study in the ID 
system. Therefore, both tripod- and A-level configura- 
tions can be used, in principle, in the realization of the 
Dirac equation. Compared with the tripod configuration 
d, 1^ , a large detuning is necessary in the A— configura- 
tion. However, the laser beams are simpler in the A-level 
configuration. Furthermore, the pseudospins in the A— 
configuration would be more robust against the collision 
of atoms since they are constructed by the lowest two 
dressed states, while the two dark states in the tripod 
configuration are not the ground states. 

We assume that the interaction can be described by 
an effective ID interacting strength g — 2h'^asN / {ml'\), 
where is the scattering length, N is the particle num- 
ber, and l± is the oscillator length associated with a har- 
monic vertical confinement. The interaction between the 
atoms (per particle) should be much smaller than the 
confinement frequency (about kHz)|20j|, and thus is also 
much smaller than Q, (in the MHz range), therefore the 
interaction can not pump the atoms outside of the near- 
degenerate subspace. Under the condition py <^ tiKy, we 
can safely neglect the Py term. In addition, we assume 
that the bosonic atoms are condensed into a BEC state. 
Within the Gross-Pitaevskii formalism, the interacting 
bosons in the near-degenerate subspace are then effec- 
tively described by a ID NLDE as ihdt'^ = Hnd'^ [ll, 
where 



HnD = -iflVyGydy 



- 7,CT^ + g^-t • + V^T + Vfc (2) 



with Vy being the effective speed of light and 7^ as the 
effective rest energy of the cold atoms. It is a remark- 
able feature that all parameters, Wy, 72 and g, can be 
controlled experimentally, providing us with a tunable 
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platform for exploration of the relativistic quantum ef- 
fects. 



III. KLEIN TUNNELING OF A SINGLE ATOM 

We now address the relativistic quantum tunneling 
that can be observed with cold atoms. To get an intu- 
itive physics picture, we first consider a single atom with 
energy E scattered by a square potential with width L 
and potential height Vs. Such a potential can be exper- 
imentally formed by a laser beam with a flat-top profile 
[l6j . The transmission coefficient Tq for the so-called 
KT regime Vg > E + "/^ , can be obtained explicitly as 

1n2 , 



Td = [l + irj-7j-ysm'{PL)/4 



where ry 



(y3--E+7.)(-E+7z) 

(£;-v;+7,)(7,-£;) 



and 



(3) 

Compared with the well-known property in nonrelativis- 
tic quantum mechanics that the transmission coefficient 
decreases mono-exponentially with the height Vs or width 
L, a distinctly different feature within this KT region is 
that the tunneling amplitude is an oscillation function of 
Vs or L even when the kinetic energy of the incident par- 
ticle is less than the height of the barrier. This relativistic 
effect can be attributed to the fact that the incident par- 
ticle in a positive energy state can propagate inside the 
barrier by occupying a negative energy state, which is 
also a plane wave aligned in energy with that of the par- 
ticle continuum outside. Matching between positive and 
negative energy states across the barrier leads to high- 
probability tunneling. We take the atoms of ^Li as an 
example. If we choose the practical parameters Ky = 10^ 
m-i, = 0.8 X 10^ m-\ fl = 10^ Hz and A = 10^ Hz, 
it is found that the Klein regime corresponds to the Rabi 
frequency fll > 0.162 MHz, which can be easily achieved 
in experiments. So we have demonstrated from a sim- 
ple example that it is feasible to observe KT with cold 
atoms. 



IV. KLEIN TUNNELING OF ATOMIC 
CONDENSATES 

As for a practical experiment it is required to release 
two conditions: the trajectory of a single atom is hard 
to detect, and it is much easier to measure the density 
evolution of an ensemble of atoms in experiments. Com- 
pared with the square potential, a Gaussian potential 
V^{y,v) — vVoe^y 1"^ , where Vq is the height and a 
characterizes the spatial variance, is much easier to be 
generated. Here v donates a barrier {v = or a po- 
tential well (y — — ), and the potential barrier (well) 
can be realized by focusing a blue- (red-) detuned far- 
off-resonant Gaussian-shaped laser beam. However, the 
conditions of resonant transmission vary with the veloc- 
ity and the width of the potential, and thus both the en- 
semble of atoms and the Gaussian potential may smooth 



the oscillations in the transmission coefficient. So it is 
natural to ask whether KT can still be observed in an 
ensemble of atoms. Surprisingly we illustrate below that 
KT of a BEC may be observed very clearly. 
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FIG. 2; (Color online), (a) A schematic diagram showing 
four kinds of scattering events, (b) Normalized density distri- 
bution in a scattering process at time t = 0, 1.5 and 2.0 ms. 
The peaks at y = are the Gaussian barriers. 

We assume that a BEC consisting of ^Li is initially 
trapped in a harmonic trap which moves along the y axis. 
At the initial time t = 0, the center of the trap is located 
at J/ = — d, and the center of the Gaussian potential is at 
?/ = 0. The trap is turned off at t = and then we calcu- 
late the evolution of the density profile of the atomic gas 
after a long enough time for scattering. The single-atom 
dispersion described in Eq. ^ is characterized by two 
branches E±{ky) = ±(7^ + h^Vyky)^^^ , where the lower 
(upper) branch represents the negative (positive) energy 
state. One can prepare an initial BEC with a designated 
mode fco at the positive or negative energy branch. The 
two branches allow us to study a more fruitful tunneling 
problem: there are four classes of scattering which de- 
scribe the wave function Vl/^ [fJ'{— ±)] scattered by the 
potential Vf{y,i'), as shown in Fig. 2(a). 

The BEC in a harmonic trap can be well described 
by a Gaussian wave packet, so we may choose the initial 
wave function as 



^/.(2/,0) 



1 



_^ifj.koy^~{y+dy/2l^o 



(4) 



where Iq is the width, fco is the central wave number 
of the wave-packet and the spinors 0^ are defined as 
0+ — (icos^, — sin^)'^, (f>- = (— i sin^, cos^)-^'' with 
^ — ■^aictan{hvyko/jz) and Tr the transposition of 
matrix. This wave function describes a Gaussian wave 
packet with the central velocity h(Ky + ^k())/m moving 
along the y-axis. After evolution governed by the Dirac- 
type Eq. ^ with time i, the finial wave function becomes 



= T exp 



HNDdt]-^^{y,0), (5) 



where T denotes the time ordering operator. We numer- 
ically calculate '^^{y^t) in Eq.® by using the standard 
split-operator method. According to the method of [l9| . 
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Eq. ([5]) can be rewritten as 

X e"K-K(2'^'')+9l*^'(2''*)l1'5* (6) 

In the sufficiently short time step St, the high-order term 
0{dt^) (due to the non-commuting) can be safely ne- 
glected. Combining with the Fourier transform between 
the position and momentum spaces, we can finally get the 
numerical solution of ^'^(y,^) following the computation 
procedure step by step with time step St. 

We have numerically calculated \E'^(?/,t), and found 
the existence of stationary solution for the scattering pro- 
cess, with an example being shown in Fig. 2(b). After 
tunneling, the incident wave packet divides into left- and 
right-traveling wave packets, and only the latter one is on 
the transmission side of the barrier. Thus we can define 
the transmission coefficient of the incident wave packet 
^ti{y,0) scattering by a potential Vf{y,i') as 

/■OO 

T^u= / *t(y,T)vI/^(y,T)d2/. (7) 

Here r (being slightly larger than cI/vq) represents a 
typical time that the reflected and transmitted wave 
packets are sufficiently away from the Gaussian poten- 
tial. One can directly measure the transmission coeffi- 
cient in Eq. ([7]) since the spatial density distribution 
Pfj^iy^'^) = [^u.(y^T)\-^ can be detected using absorption 
imaging |17| . 

We first look into the tunneling phenomena for a BEC 
in the absence of interactions {g = 0). We note that there 

are two identities T_|_+ = T__ and T y ~ since 

Eq.@ with g = is invariant under the charge conjuga- 
tion ^18|]. We plot the transmission coefficient T++ as a 
function of the height Vq and width a in Fig. 3(a) with 
the practical parameters. It is interesting to note that the 
transmission coefficient decreases exponentially to zero 
with Vg when Vq < Vq , while it increases and then 
is an oscillating function in the Klein region Vq > Vq , 
these results are similar to the results of Eq. for the 
square barrier. Here the critical value of the potential 
height may be estimated approximately using the square 
barrier with = E{ko) + ^ 0.09 MHz. Moreover, 
the feature T++ = T__ is also confirmed in the inset in 
Fig. 3(a). As for the transmission coefficient T+_|_((t), we 
may obtain several tunneling oscillations with the poten- 
tial width, but it decreases to when the width is further 
increased. Although the amplitude of tunneling oscilla- 
tion is less than the unit compared with the tunneling of 
a single atom, the amplitude of tunneling oscillation can 
be more than 0.5 and meanwhile the period can be a few 
micrometers, which is experimentally detectable. 

Another interesting feature induced by relativistic ef- 
fects is that, a BEC with negative energy can almost 
completely transmit a wide Gaussian potential barrier, 
as shown in Fig. 3(b). The transmission coefficient T y 



is an oscillating function of the potential width a when 
a is smaller than 3 fim, while it saturates quickly to the 
unit when the potential width is larger than 3 fim, lead- 
ing to the unexpected result that a wide Gaussian po- 
tential barrier is actually totally transparent for a BEC. 
This phenomenon can be understood through the fact 
that this scattering feature is actually equivalent to that 
of a BEC of positive energy scattered by a Gaussian po- 
tential well because of T_+ = T+_. We also calculate 
the transmission coefficient for the central mode of the 
wave packet, as shown in the inset in Fig. 3(b), which 
further confirms that a wide enough Gaussian potential 
well is transparent. The reason lies in the fact that, in 
contrast to the periodic function (without a saturation 
value) in a square potential well, the Gaussian potential 
well is smooth in the whole space, and thus can even sup- 
port adiabatic motions of wave packets in the large width 
limit. 




FIG. 3: (Color online). KT of BECs. (a) T++{a) for 
Va/h = 0.2 MHz and T++{Vg) (inset) for cr = 5 ^m. The 
tunnelings of a BEC with the classic kinetic energy term and 
the conventional atomic interaction (A'^ = 2x lO'', l± = lA nm 
and Us = 5ao with ao being the Bohr radius) are also depicted, 
(b) Coefficients T ^{a), {T{a)) (insert) of one atom with cen- 
tral mode ko, and of 10"* atoms for Vo/ti = 0.2 MHz. The 
other parameters in (a) and (b) are lo = 10 fim, ko = 5.5 x 10^ 
m"\ 7^/^ = 30 kHz, and d = 4{lo + a). 




a (nm) a (urn) 



FIG. 4: (Color online). Comparison of KT of BECs with 
that of an ensemble of non-condensed atoms, (a) T++{a) with 
iV = 10^ 4 X 10^ 10" are shown by fixing the energy Eint for 
lo = 10 fim and d = 4(;o-|-o-). (b) (T(cr)) with Na = 1, 10, 10^, 
and 10'* atoms for at = 5 x 10^ m~^. The other parameters 
in (a) and (b) are Vc/h = 0.2 MHz, fco = 5.5 x lO'^ m"\ and 
= 30 kHz. 



5 



The tunneling properties exhibited in Fig. 3(a) and 
3(b) are intrinsic relativistic and macroscopic quantum 
phenomena that can not be explained with an incoherent 
ensemble average of many atoms. To clarify this point, 
we calculate the average transmission coefficients for an 
ensemble of Na noninteracting atoms defined as 

(T) = (8) 

^ i—l 

where T(ki) donates the transmission coefficient for atom 
i with the wave number ki scattered by the potential. 
The numerical calculation method for T(ki) is given in 
the Appendix. Here we choose ki to be the same Gaus- 
sian distribution as that of the initial BEG wave function 
^Ai(y,0), i.e., ki ^ N{ko,al) with the variance Cfc = l/lo- 
The (T) of 10^ atoms is shown in the inset in Fig. 3(b), 
which is almost the same as that of a single atom since 

(Tfe is small. The differences between (T) and T [_ in 

Fig. 3(b) demonstrate that the tunneling of BEG is not 
equivalent to an ensemble average of the individual atoms 
even with the same distribution of wave number. The co- 
efficient (T) represents an incoherent transmission of the 
individual particles since it is a sum of the transmission 
coefficients of all particles. In contrast, the phases of all 
atoms in the BEG are the same and then the transmis- 
sion of a BEG is coherent. The coherent transmission in 
a BEG and incoherence in (T) cause the difference in Fig. 
3. All particles in the BEG are in the same phase because 
the macroscopic number of particles are condensed in the 
same state, so the coherent transmission of the BEG may 
be called macroscopic quantum tunneling. 

To clarify further the macroscopic quantum phenom- 
ena in the relativistic tunneling of a BEG, we compare 
the scaling properties of transmission coefficients for a 
weakly interacting BEG and an incoherent ensemble av- 
erage of atoms. An example of scaling of T++ is plotted 
in Fig. 4(a). In the calculations, we have fixed the weak 
interatomic interaction energy Ei„t ~ g/lo and kept the 
parameter 7 = mglo/NH^ 1 [20] for Iq = 5 /im when 
N = 10^, l± = 1.4 ^m and = 5ao (7 ~ 10"^), both 
of which restrict our discussions in the regime for ID 
BEGs, where Dirac dynamics instead of nonlinear dy- 
namics dominates. In this case, the increase in parti- 
cle number is achieved by proportionally increasing the 
length /q of the BEG with small 7. For comparison, we 
also calculate the scaling of (T) for the atom numbers of 1 
(with ki = fco), 10, 10^, and 10^ in Fig. 4(b). Gomparing 
Fig. 4(a) with Fig. 4(b), a distinct difference between the 
BEG and the ensemble average of the individual atoms is 
that, the coefficient r-|_+ increases with increasing atomic 
number of the BEG, while the coefficient (T) decreases 
with the increasing of the atomic number. However, in 
order to keep the same interaction parameter in the above 
calculation, we have increased simultaneously the parti- 
cle number and the width of the Gaussian wave packet. 
In this way the momentum distribution of the wave func- 
tions is shrunk, which is a dominant reason for the above 



scaling feature. That many atoms may condense into the 
same momentum state is essential for the observation of 
KT in a BEG. 



V. DISCUSSION AND CONCLUSION 

Before concluding, we wish to make two additional 
comments, (i) To judge the feasibility of the Dirac ap- 
proximation in Eq. ^ , the coefficients T_|_+ with or with- 
out the quadratic term are compared in Fig. 3(a). It is 
shown that the quadratic term leads to merely a slight 
left-shift of the tunneling peaks. This phenomenon can 
be interpreted by the fact that the wavelength of the 
BEG inside the barrier decreases slightly in the presence 
of the additional low kinetic energy. This result verifies 
that the approximation leading to the Dirac equation 
is well satisfied, (ii) In Fig. 3(a), we have also calcu- 
lated the transmission coefficient for BEGs with conven- 
tional atomic interactions without Feshbach resonance, 
in which case the experimental setup can be simplified. 
The result shows that the effect of the realistically weak 
interaction is small; it merely smooths the tunneling os- 
cillation slightly. Therefore the exotic tunneling phenom- 
ena addressed here survive in the case of weak interaction 
between atoms. 

In summary, we have proposed an experimental scheme 
to detect macroscopic KT using a spin-orbit-coupled 
BEG. Through numerical simulations, we have elabo- 
rated that such macroscopic KT can be observed under 
realistic conditions. In view of the fact that a spin-orbit- 
coupled BEG was realized in a very recent experiment 
[l3|], it is anticipated that the present proposal will be 
tested in an experiment in the near-future. 
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APPENDIX: The derivation of T(fc,) in Eq.dH) 

The Dirac equation for particle scattering by a Gaus- 
sian potential can not be solved analytically for an in- 
coming atom with energy Ei = ^ (hvyki)"^ + 7^ and mo- 
mentum Pi — hki. However, here we adopt an effi- 
cient method to solve it numerically based on transfer 
matrix methods Q. The numerical procedures are out- 
lined as follows. First, one cuts the Gaussian potential 
into a spatially finite range y £ [— j/c,2/c], where the cut- 
off position j/c should be chosen to guarantee that the 
potential height outside the range is low enough to be 
transparent for the atoms, i.e., V^{yc) ^ Ei^Vo- Sec- 
ond, one divides this range equally into n spindly seg- 
ments, and each segment may be considered a square 
potential if n is large enough. The potential height of 
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the j-th (j = 1,2, •••,«) square potential is given by 
= y^iVj + //2) with = -y, + (j - 1)/ and the 
width of each potential / = 2yc/n. In this case, the 
Gaussian potential can be viewed approximately as a 
sequence of connective small square potential barriers, 
and thus the transmission coefficient T{ki) k, l/jmnp, 
where mn is the first element in the whole transfer ma- 
trix M = MnMn-i ■ ■ ■ Mj ■ ■ ■ M2M1. Here Mj denotes 
the transfer matrix of the j-th square potential barrier, 
whose explicit elements are given by |9| 

(M,)n = (cos£^+z^±;isin£|i)e-iP-/, 

(^.)i2 = (*4i^sin£^)e-^'-(^^+^^+^), (9) 
(M,)2i = {M,)l2, 

(M,)22 = (M,)ti, 



where k = {Ei-jz)/{vyPi) and Kj = {Ei-^^-Vj) / {vyPj) 
with [Ei — VjY — VyP^^ + 7^. Note that this numerical 
calculation scheme recovers the non-relativistic scatter- 
ing governed by the Schrodinger equation. 
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